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Context. Whenever stars are rotating very fast (fi/n cr it > 0.7, with fi C rit the Keplerian angular velocity of the star 
accounting for its deformation) radiative stellar winds are enhanced in polar regions. This theoretical prediction is now 
confirmed by interferometric observations of fast rotating stars. 

Aims. Polar winds remove less angular momentum than spherical winds and thus allow the star to keep more angular 
momentum. We quantitatively assess the importance of this effect. 

Methods. First we use a semi-analytical approach to estimate the variation of the angular momentum loss when the 
rotation parameter increases. Then we compute complete 9 Mq stellar models at very high angular velocities (starting 
on the ZAMS with fl/fi cr it = 0.8 and reaching the critical velocity during the Main Sequence) with and without 
radiative wind anisotropics. 

Results. When wind anisotropics are accounted for, the angular momentum loss rate is reduced by less than 4% for 
Sl/Q CI it < 0.9 with respect to the case of spherical winds. The reduction amounts to at most 30% when the star is 
rotating near the critical velocity. These values result from two counteracting effects: on the one hand polar winds 
reduce the loss of angular momentum, on the other hand, surface deformations imply that the mass which is lost at 
high co-latitude is lost at a larger distance from the rotational axis and thus removes more angular momentum. 
Conclusions. In contrast with previous studies, which neglected surface deformations, we show that the radiative wind 
anisotropies have a relatively modest effect on the evolution of the angular momentum content of fast rotating stars. 
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1. Introduction 

Over the last years, we added to the Geneva stellar evolu- 
tion code several major improvements, as rotation, the in- 
clusion of extended nuclear reaction network allowing to ex- 
plore the advanced phases of massive star evoluti on (neon, 
oxygen and silicon burning, see lHirschi et al.ll2004r) . and the 
inclusion of magnetic field . 

The inclusion of rotation improved the agreement be- 
tween the outputs of numerical models and observational 
results, as the surface enrichments, the ratio of blue to red 
supergiants in the SMC, the ratio of WR to O-type stars 
and the ratio of type Ibc to type II supernovae (see, for 
example iMevnet fc Maederl 120001: IMaeder fc Mevnetll200ll 



IMevnet fc Maederl l2003t iGeorgy et al.l I2009D . The treat- 
ment of internal magnetic field in the models gives a more 
realistic rat e of Gamma Ray Bursts as a function of the 
metallicity (|Yoon et al.ll2006t ). produces a rotation rate for 
the you ng pulsars in much better agreement with the obser- 
vations ( Heger et al.ll2005l ), and allows to explain t he flat 
rotation profile of the Sun (Eggcnb crger et al.ll2005[ ). 

Rotation acts not only in the interior of the stars, 
however, but has also several effects on the surface. 
IMaeder fc Mevnetl (|2000f ) show that rotation increases the 
global mass loss rate. It also modifies the shape of the star, 
and consequently surface quantities such as: effective grav- 
ity, effec t ive te mperature, and radiative flux. Interestingly, 
iMaederl (|2002t l shows that the mass loss in fast rotating 
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massive stars does not remain isotropic, but becomes in- 
creasingly anisotropic as the rotation approaches the criti- 
cal rotation parameter. They show that this favours a bipo- 
lar stellar wind, and modifies the quantity of angular mo- 
mentum removed from the star. 

With the developments of the interferometric technics, 
some of the predic ted effects have recen tly become observ- 
able. For example, ICarciofi et al.l (|2008j ) have obtained for 
the ratio of the equatorial to polar radius a value of 1.5 for 
the very fast rotating star Acher nar as it is expected i n the 
frame of the Roche model theory. iMonnier et al.l (|2007t ) pro- 
vide a map of the effective temperature over the surface of 
Altair, showing that the temperature gradient between the 
pole and the equator of t his star is in good agreement with 
the vo n Zeipel theorem (jvon Zeipell 119241) . iMeilland et all 
(|2007t ) observed evidences of a disc and of a polar wind 
around the star /i Ara, as is predicted for a star rotating at 
the critical limit. 

At the moment the effects of wind anisotropies on stel- 
lar models have been quantitatively explored only in two 
previous publications (jMevnet fc Maederl [2001 12007( 1 . In 
the present work we improve the numerical treatment over 
these works and reexamine this effect. We include the fol- 
lowing improvements: 

— For the first time, we account not only for the variation 
of the mass flux with the co-latitude as was done in 
the previous work, but we also account for the surface 
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deformation of the star. As we will see this last effect 
cannot be neglected. 

— W e use an update d expression for the mass flux obtained 
bv lMaederl (120091) . 

— We accounted for the variation of the force multiplier 
parameters over the surface of fast rotating stars (see 
below for more details on that point). 

To check the validity of the numerical results and also 
to study in a clearer and simpler context the various effects 
intervening in the loss of angular momentum, we propose a 
semi-analytical approach to estimate the effect of fast rota- 
tion on the loss of angular momentum. The relative effects 
obtained in that way depend only on one parameter, the 
ratio Q/f2 cr it, where fl is the surface angular velocity and 
f2 C rit) the critical angular velocity, i.e., the angular velocity 
such that the centrifugal acceleration at the equator com- 
pensates for the gravity (accounting for the deformation of 
the shape of the star). 

The paper is organised in the following way: in Sect. [5J 
we give the theoretical aspects of wind anisotropy. The 
third section presents our semi-analytical approach. In 
Sect.Sl we discuss results based on complete numerical stel- 
lar models. Conclusions are presented in Sect. [5j 



2. Rotation and wind anisotropy 

2.1. Increase of the global mass loss rate induced by rotation 

As shown by|Maedcr fc Mevnetl (|2000h . the local radiative 
mass loss rate AM by unit surface Act can be written: 
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In this expression A = (ka)<* f^- 2 ) a , where a and fc 
arc the force multiplier parameters empirically determined 

(jLamers et al.lll995l) . M* = M (l - 2 ^ p j is the reduced 

mass, where p m is the internal average density, L is the stel- 
lar luminosity, g e g is the effective gravity at the co-latitude 
9 (i.e., the vectorial sum of the gravitational acceleration 
and of the centrifugal one) and ((9) expresses the deviation 
from the von Z eipel theorem produced by shcllular rotation 
(|Maederlll999l this term is generally negligible), Tq(9) is 
the local Eddington factor, taking into account the effect 
of rotation: 
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with k cs the electron scattering opacit y. The term 1/ 8 in 
the power of expression ((T|) was added b vlMaederl (120091 see 
his chapter 14.4) and does not appear in lMaeder fc Mevnetl 

(|2000L a T cS 2 was absent in their expression 4.24). Note that 
in the frame of the line driven wind theory used here to ob- 
tain the expression of the mass flux (Eq. (TTJ)), the opacity is 
exp ressed in as a funct ion of the electron scattering opac- 
ity (jCastor et al.|[l975l ). The variations are accounted for 
through the force multiplier parameters a and fc. 



Averaging expression (TTJ) over the whole stellar surface 
E, one obtains the total mass loss rate of the star 
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This allows us to compute the ratio of the mass loss rate 
of a rotating star to the mass loss rate of a non-rotating 
one lying at the same position in the Hertzsprung-Russel 
diagram: 
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where Lsdd is the classical Eddington factor for a non- 
rotating star. As a result, we see that the faster the star 
rotates, the more mass will be lost per time unit. 

2.2. Wind anisotropy 

According to lMaede"n (|1999f ). two main effects contribute to 
the development of anisotropies in the stellar winds. The 
first, called the g e ff-effect, is due to the variation of the 
effective gravity with the co-latitude: g c g is smaller at the 
equator than at the poles, and thus, the mass loss, which is 
directly related to g e g (see eq. (HJ ) , is favoured at the poles 
for a rotating star. 

The second effect is called the re-effect. Due to the so- 
called bistability in the stellar winds (see lLamers et al.l 
119951) . the A term in eq. (UJ increases for lower values of 
the effective temp erature, i.e., towards the equatorial re- 
gions (see Fig. 6 in lEkstrom et al.ll2008|) . This occurs when 
the effective temperature are below 11500 K (Lamers, pri- 
vate com.). This favours an equatorial mass loss. 

Looking at cq. ([T]), wc also expect a contribution to the 
latitudinal variation of the mass loss due to the term n{9) 
in r n (0). 

2.3. Critical velocities 

As the concept of critical velocity is treated in very dif- 
ferent ways through the literat ure, we briefly repe a t som e 
common definitions. Following iMaeder fc Mevnetl (|2000D . 
we define two critical velocities. The first one is the tra- 
ditional Keplerian velocity at the equator when the star 
rotates at the critical velocity defined by g c g = 0: 
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where i? e b (Rpb) is the equatorial (polar) radius when the 

first critical velocity is reached. The numerical factors 2/3 

comes from the polar to equatorial radius ratio when the 

star is at the criti cal velocity and the R oche approximation 

is valid (see, e.g.. lEkstrom et al"1l2008l ). We also define the 

critical angular velocity f2 cr it = -fp^ , and the ratio ui = 

n cb 
n crit • 

The second critical velocity is reached when the star is 
at the so-called fir-limit, i.e. when the local Eddington fac- 
tor (accounting for the effects of rotation) defined in cq. © 
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is 1. According to iMaeder fc Mevnetl (|200C 
equal to: 

81 1 - r Edd GM d2 

-rt„, 



v 2 
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this term is 
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with TEdd the Eddington factor, Vb is the ratio of the vol- 
ume enclosed by the surface when the star rotates at the 
second critical velocity to the volume of a sphere with a 
radius equal to R p b , R c is the actual equatorial radius. 
For r cd d < 0.639, u crit ,2 = «crit,i- For r Ed d > 0.639, 
^crit,2 < ^crit.i, arid the relevant critical velocity is the sec- 
ond one, because it is reached first. 



3. Effects of fast rotation on angular momentum 
loss: a semi-analytical approach 

In this section we derive the variations of 

— the shape of the star, 

— the radiative mass flux, 

— the angular momentum flux, 

— the global angular momentum loss, 

as a function of only one parameter: uj. Said in other words, 
the results obtained are independent of the mass, metallic- 
ity, and evolutionary stage of the star considered. To obtain 
such a simple dependence, some hypothesis and normaliza- 
tions have to be made: 

— We use the Roche approximation for computing 
the gravitational p otential (which is valid here, see 
IMevnet et aLI l20lj) . Note that we also assume the 
Roche approximation in our numerical stellar models. 
This assumption is not only supported by the direct 
observations of (a few) rapidly rotating stars, but also 
by the fact that, in our numerical models, the angular 
velocities are far below the critical velocity in a large 
fraction of the total mass. The results presented in this 
work are however probably dependent on the validity of 
this hypothesis. 

— We neglect the variations over the surface of the force 
multiplier parameters (see below). This variation will be 
accounted for in more complex stellar evolution models 
(see Section [43)) . 

— We neglect the correcting factor ((8) in Eq. [TJ 

— The angular velocity of the surface is assumed to not 
depend on the colatitude 8 (no differential rotation of 
the surface). 

— We suppose that the polar radius of the star remains 
constant when uj increases from (no rotation) to 1 
(critical rotatio n). This is well verifi ed in complete nu- 
merical models (jEkstrom et al.ll2008l ). We normalise the 
polar radius R p to 1. 

— The total mass loss rate is taken equal to An for all 
rotation velocities. In that case, the mass loss rate per 
unit surface is equal to 1 in the non-rotating case. 

3.1. Shape of the surface 

With the above hypothesis, the equation of the surface of 
the star can b e given as a function of the rotation rate uj 
(lMaederll2002h 




Fig. 1. Shape of the surface for various values of uj = j^— 
(labelled at the bottom of each curve). The x-axis is the 
equatorial radius, and the y-axis the polar one. Hence this 
is how we would see the star equator-on. 



where x{uj, 8) 



r(u>,8) 
Rr> 



is the ratio of the radius at a given 
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co-latitude to the polar one. We can easily express 8 as a 
function of the normalised radius x: 



if uj ^ 
ifw = 0. (8) 



The range of satisfactory values for x is a function of uj. 
It starts from 1 (to have a positive value under the square 
root), and goes up to the first positive root of the equation 
4uj 2 x 3 — Tlx + 27 = 0. In Fig. [TJ we show how the shape 
of the surface varies for various values of w, starting from 
oj = (non rotating case) to uj = 1 (critically rotating case) . 
As uj increases, the centrifugal force deforms ever more the 
star, and the equatorial radius increases. When the star is 
exactly at the critical angular velocity {uj = 1), we see from 
eq. ([JJ that the equatorial radius is 1.5 times larger than 
the polar one. 

Eq. [71 which allows to deduce the shape of the star, 
depends only on uj. Maximal deformation is obtained for 
<jj = 1. Stars reaching the second critical velocity will have 
a ratio uj below 1. Despite being at the critical velocity, they 
will not show as strong deformations as stars reaching the 
first classical critical limit. Since, as we will see, it is the de- 
formation of the star which triggers the wind anisotropics. 
This means that stars which would be at the fir-limit do 
not present as strong wind anisotropies as stars at the clas- 
sical fi-limit. Accordingly, the fact that 77-Carinae presents 
strong polar winds, implies that this star should rotate at 
velocities close to the first classical critical rotation veloc- 
ity, in case the bipolar shape is due to rotationally induced 
wind anisotropy. 
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Fig. 2. Left panel: Mass loss rate per surface unit as a function of the co-latitude 9. The pole is on the left, the equator on 
the right. The mass loss rate per surface unit is normalised to 1 in the non-rotating case, and thus is expressed without 
unit. Right panel: 2D colour plot of the mass loss rate per surface unit. The x-axis is the co-latitude, the y-axis the 
rotation parameter uj — J7/J7 cr ; t . The colour scale on the right indicates the local mass loss rate per surface unit: from 
blue for the lower mass loss flux to red for the higher. 



3.2. Mass flux variations with the latitude 



with 



To study how the mass flux is modified by the rotation, 
we use eq. [T] to compute the local mass loss rate per unit 
surface. As said above, we normalise the total mass loss to 
a value of 47r, in order to have a local mass flux of 1 at 
every co-latitude in the non-rotating case. 

The results of these calculations are shown in Fig. [5] 
Not surprisingly, the mass loss rate per surface unit is con- 
stant when there is no rotation. As the rotation param- 
eter increases, the mass flux increases towards the pole, 
and decreases towards the equator, due to the variation of 
the effective gravity, producing a strong anisotropy in the 
winds. Typically the mass flux at the pole is greater than 
two times the mass flux at the equator for u> > 0.8. In the 
extreme case, when the rotation parameter u = 1, the effec- 
tive gravity at the equator is zero, and the radiative mass 
flux becomes also zero in this idealised representation. 

From Fig. [2j we see that the faster the star rotates, the 
more mass will be lost in the polar region. Integrating the 
mass flux from the pole to a given co-latitude 9 and dividing 
by the total mass loss rate, we obtain the fraction of the 
total mass flux in a cone of semi-aperture 9, for a given 
rotation factor w: 



M . 



(9) 



where da (8) is the surface clement at the colatitude 8, given 
by da = r ^^JfJ^^ \ with e the angle between the local 
effective gravity and the radial direction. The angle e is 
computed using the components of the effective gravity g e g: 

9eff • e r 



cos(e) 



I9cff| 



(10) 



+ n 2 r(9)sm(9)cos(9)e g , (11) 

where e r (eg) is the radial (colatitudinal) unit vector. 

The result is shown in Fig. [3] For a non-rotating star 
(lower curve), we see that 50% of the total mass is lost in 
a cone of semi-aperture 60°. When the star is at the first 
critical velocity, the aperture of the cone containing 50% of 
the total mass flux is slightly reduced: it is around 48°. 

3.3. Latitude dependency of the angular momentum loss 

Once we know the local mass flux and the shape of the 
surface, it is possible to compute the local loss of angular 
momentum induced by the stellar winds, for a given angular 
velocity of the surface (let us recall that we suppose in this 
work that the angular velocity of the surface Q, is constant 
over the whole stellar surface). 

The loss of angular momentum per surface unit and 
time is given by: 



dC 
da 



AM 
~Aa~ 



(9)nK 2 (9), 



(12) 



where 1Z(9) is the distance from the considered unit sur- 
face element to the rotation axis at the co-latitude 9. Using 
the surface element (see above), and integrating over ip to 
obtain only the co-latitudinal variation of the angular mo- 
mentum loss, we have: 



dC n AM 
d9= 2n ^a- 



(9)r 4 (9)n 



sm 3 (9) 
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20 40 60 80 

B[°] 

Fig. 3. Fraction of the total mass loss contained in a cone 
with a semi-aperture 9 given by the x-axis. The lower 
curve is the result for a non-rotating star losing its mass 
isotropically. The upper curve is for a critically rotating 
star (lj = 1). The intermediate dashed curves are the re- 
sults of the same rotation parameter omega as in Fig. [2j 
i.e. lo = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, 0.995 respectively 
(from bottom to top). 



which is the contribution to the total angular momentum 
loss of an infinitesimal ring at the co-latitude 9 and extend- 
ing over an angle d9. In order to avoid the fi-dependency, 
and as in our model, the surface of the star rotates at a con- 
stant angular velocity, we consider further the ratio C/fl. 
This permits to easily compare models with various rota- 
tion parameters. The distribution of the ratio C/Q brought 
away by the wind is shown in Fig. 0] 

In this figure, we see, how the angular momentum flux 
(normalised by the surface angular velocity) is distributed 
as a function of the co-latitude 9 for various rotation param- 
eter oj. Two effects are in concurrence: first, the increase of 
the equatorial radius (see FigO}, which increases the angu- 
lar momentum flux near the equatorial regions, and second, 
the decrease of the local mass loss rate near the equator, 
which decreases the angular momentum flux in the same 
area. 

Without rotation, the mass loss rate per surface unit 
is constant over the whole surface of the star, and the 
angle e between the effective gravity direction and the 
radial direction is zero. Examining eq. (|13[) , we see that 

d {j^j /d# varies as sin 3 (9) (since all other terms are con- 
stant). The corresponding curve is labelled ui = in Fig. [4] 
Progressively increasing the rotation parameter, we see that 
the deformation of the stellar surface produces an increase 
of the angular momentum loss in the equatorial region. 
Once the rotation parameter u> ~ 0.75, the increase of the 
equatorial radius becomes counterbalanced by the progres- 



sive decrease of the local mass loss flux in the same region. 
The angular momentum loss becomes thus more and more 
reduced at the equator, and the maximum of momentum 
loss is shifted towards the pole, up to a co-latitude of ~ 70° 
when the rotation becomes critical. At that moment, no 
more angular momentum is lost at the equator, since there 
is no equatorial mass loss through radiative winds in this 
regime in our model. 

3.4. Total angular momentum lost 

Once the distribution of the angular momentum loss is 
known, we can integrate it over the co-latitude to obtain 
the total angular momentum loss rate due to the stellar 
winds: 

Jo Aa u w cos(e) v ; 

To well understand the effect of the anisotropic winds on 
the total angular momentum loss, we distinguish the fol- 
lowing cases: 

— Case 1: we determine a mean stellar radius r mcan using 
the following relation: 

L = £ CT <T e 4 ff ) ee Anr 2 mean c7 (T c %) (15) 

where L is the stellar luminosity and S the total stellar 
surface. We neglect the stellar deformation and the wind 
anisotropy, and the loss of angular momentum is thus 
computed on a sphere of radius r moan : £ = |MQr^ eall . 
This radius is the one we would find if we measure the 
luminosity and the effective temperature of the star, and 
suppose that it is perfectly spherical. It gives the angular 
momentum loss as computed in numerical models where 
the effects of rotation on the shape of the surface and 
the mass loss distribution is neglected; 

— Case 2: the deformation of the star is accounted for, 
but the mass loss is uniformly distributed over the stel- 
lar surface (i.e. the winds anisotropy is not taken into 
account). This case is academical, but is interesting in 
the sense that it allows to see the effect of the deforma- 
tion only; 

— Case 3: both the deformation of the shape of the sur- 
face and the anisotropy are accounted for. The angular 
momentum loss is computed with relation [T4l 

The results we obtain are shown in Fig.[SJ The top panel 
shows the variation of the ratio C/VL removed by the wind in 
the cases "neither deformation, nor anisotropy" , computed 
with the mean radius discussed above (£/f2 sp h, case 1), 
"deformation only" (£/f2i so , case 2), and "deformation + 
wind anisotropy" (£/fi an i, case 3). In the three cases, the 
value for u> = is given by the integration of eq. [Ml with 
AA//A<7 = 1, r = 1 and e = at each colatitude, leading 
to C/fl = 271-/* sin 3 (6»)d6» = 8tt/3. 

The increase of £/il sp h in case 1 is entirely due to the 
increase of the mean radius r moan when u> increases. When 
the deformation is accounted for, we see that more angu- 
lar momentum is lost because most of the mass leaves the 
surface of the star at a greater distance from the rotational 
axis. We see that, at the critical limit, deformation would 
increase the angular momentum loss rate by around 20% 
with respect to case 1. For case 3, we see that the wind 
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Fig. 4. Left panel: Rate of angular momentum loss (normalised by the surface angular velocity) from infinitesi- 
mal rings centred on the rotational axis for each co-latitude 9. The width of each ring is given by rdO. The ro- 
tation parameter u> is indicated for the curves with uj = and u = 1. The intermediate dashed curves are for 
io = 0.5, 0.75, 0.9, 0.96, 0.98, 0.99, 0.995 respectively. The rate of angular momentum loss is expressed in terms of the 
normalised units M+R^/ (Air) where M* is the mass loss rate of the star, and i?* the polar radius of the star. Right panel: 
Same as left panel in a 2D colour plot. The cc-axis is the co-latitude, and the y-axis the ratio u> = Q/Q CI n. The colour 
scale on the right indicates the local angular momentum loss from infinitesimal rings centred on the rotational axis for 
each co-latitude 9. The units are the same as in left panel. The smaller values are in blue and the larger ones in red. 



anisotropy largely compensates for the effect of the defor- 
mation and decreases the rate of angular momentum loss 
by 49% with respect to case 2 and by 25% with respect to 
case 1. 

The lower panel of Fig. [S] shows the ratio £ sp h/£ an i. It 
shows the real impact of the account for the anisotropy of 
winds and the deformation of the stellar shape, compared 
with a model where we consider an isotropic spherical wind 
on the surface, with a radius determined by the stellar lu- 
minosity and mean effective temperature (as most of the 
stellar evolution codes). 

Interestingly enough, the error committed on the an- 
gular momentum loss when neglecting the effects of wind 
anisotropies is small in most of the cases. It is less than 
4% if 10 < 0.9. At the critical velocity, the error is larger. 
Up to 25% more angular momentum can be kept in the 
star when the effects of wind anisotropies are accounted 
for. Therefore, the effects of the wind anisotropics become 
important only for the faster rotators. This indicates that 
for most of the cases studied in stellar evolution, the precise 
account for the anisotropies are not relevant, and the errors 
induced by neglecting it will remain small. 



4. Effects of fast rotation on angular momentum 
loss: a numerical approach 

The analytic relations above can provide some orders of 
magnitude estimates of the impact of the wind anisotropies 
on the loss of angular momentum assuming that to remains 
more or less constant as a function of time. To obtain more 



accurate values, it is necessary to compute numerical stellar 
models. 

An additional complication, let aside in the above es- 
timate, comes from the fact that the global mass loss rate 
increases with faster rotation. We first consider a case where 
the mass loss is kept constant. This will allow to make a 
more direct comparison with the semi-analytical results ob- 
tained above and thus to check that the implementation of 
the process was done correctly in the stellar evolution code. 
In a second step (see Sect. 14. 3p . we shall consider the case 
of a model with all the usual prescriptions, in particular ac- 
counting for the evolution of the mass loss rate as a function 
of time, and the possible variation of the force multiplier 
parameters a and k over the stellar surface. 

4.1. Models with constant mass loss rate 

We examine the effects of the anisotropic stellar winds on 
the evolution towards the critical velocity of two sets of 
models of 9 M Q star. To see the effects of anisotropic winds 
already on the zero age main sequence (ZAMS), we start 
the computation of the stellar models at a very high ini- 
tial angular velocity of 80% of the critical velocity. The 
metallicity is taken equal to Z = 0.002, i.e., equivalent to 
that of the Small Magellanic Cloud (this implies smaller 
mass loss rates than at higher metallicity, favouring the 
reaching of the critical velocity for the numerical mod- 
els presented in S e ct. | 4.3p . The rotation is treated as in 
iMaeder &; Mevnetl (|2005l ). accounting for the internal mag- 
netic fiel d and its imp act on the transport of angular mo- 
mentum (lSpruitll2002r) . These numerical models are based 
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Fig. 5. Top panel: Evolution of the total angular momen- 
tum loss as a function of the rotation rate. The long-dashed 
curve is the angular momentum loss when the mass loss is 
supposed isotropic and spherically symmetric (C sp h) , case 
1). The solid curve is the angular momentum loss when the 
deformation of the star is accounted for, but the mass loss is 
assumed to be isotropic (Asoj case 2). The dotted curve is 
the angular momentum loss when both the deformation and 
the anisotropy of the winds are accounted for (C an i, case 
3) (see text for more details). Bottom panel: Ci so /C an i as 
a function of u. 



on the shellular rotation assumption (jZahnl [19921. Even if 
the surface of our models reach rotation parameters close to 
1 , as a large part of the stellar interior rotates far from the 
critical velocity, we consider that this assumption is valid. 
The account for the magnetic field ensures a strong cou- 
pling between the centre and the surface of the star, and 
leads to higher surface velocity. We expect thus a more im- 
portant effect of the anisotropic winds in this context. Both 
models were computed using a constant mass loss rate of 
10 -9 M Q yr _1 , independent of the stellar surface parame- 
ters, and independent of the rotation rate. One model was 
computed with the account for the anisotropic winds, and 
one assuming isotropic winds. In both cases, the deforma- 
tion of the shape of the star is accounted for. 

In Fig. [6l we show the Hertzsprung-Russcl diagram 
(HRD) of both models. The evolution stages where the 
models reach 80%, 90% and 95% of the critical angular 
velocity are indicated on the tracks. The ZAMS is bottom- 
left, and the evolution proceeds towards the top-right cor- 
ner. We see that the account for the anisotropic winds has 
only a minor impact on the evolutionary track. A small de- 
viation begins to appear when the surface velocity is around 
90% of the critical velocity. The anisotropic model evolves 
slightly more on the red side of the HRD. This is because 
this model rotates faster than the isotropic model (see the 
top panel of Fig. [7]). Its surface is slightly larger, and thus, 



t 1 1 1 1 1 1 1 1 r 




4.42 4.4 4.38 

l°g(T eff ) 

Fig. 6. HR diagram for the isotropic model (black solid 
line) and anisotropic model (red dashed line) with constant 
mass loss rate. The rotation parameter is indicated for some 
points along the tracks. The end point of the track (on 
the right) corresponds to the reaching of the first critical 
velocity. 



for a given luminosity, the mean effective temperature will 
be lower. 

Fig. [7] shows on the top panel the evolution of the ratio 
f2/O cr it as a function of time for both models. The bottom 
panel shows the total angular momentum of the star. As 
expected, the mass loss due to the stellar winds causes a 
decrease of the angular momentum kept in the star. The 
model with the account for the wind anisotropy loses less 
angular momentum than the isotropic model. As a conse- 
quence, the stellar surface of the anisotropic model rotates 
slightly faster than the isotropic model. 

The mean f2/f2 C rit fc> r the anisotropic model is f2/i7 cr jt = 
0.867. According to Fig. [SJ we expect that the anisotropic 
model keeps 1.12 times more angular momentum than the 
isotropic one. During the time between the ZAMS and 
the reaching of the critical rotation rate, the anisotropic 
star loses an amount of angular momentum A£ an ; = 
3.56 • 10 50 gcm 2 s -1 . During the same time, the isotropic 
model loses A£ iso = 4.13 • 10 50 gcm 2 s- 1 . The final ratio 
A£i so /A£ ani = 1.16 is very close to the estimate based on 
the mean f2/f2 cr it. 

4.2. Total mass loss rate as a function of 'us and T^d 

Before discussing the results of models accounting for time 
dependent mass loss rates, let us briefly recall how rotation 
enhances the global mass loss rate. Examining eq. (j4|), we 
see that the mass loss rate of a rotating star is simply ex- 
pressed as a function of the angular velocity f2, the classical 
Eddington factor T^dd and the local Eddington factor To- 
Let us introduce in this relation the rotation parameter u> 
as defined above, and the definition of given by relation 
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Fig. 7. Top panel: Rotation parameter as a function of the 
time for the isotropic (black solid line) and anisotropic 
model (red dashed line) with constant mass loss rate. 
Bottom panel: Total angular momentum contained in the 
star in unit of 10 53 gcm 2 s _1 . The models are represented 
as in top panel. 
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Fig. 8. Logarithm of the ratio M(lo)/M(uj = 0) as a func- 
tion of uj = -gp— for various values of the Eddington factor 
(value indicated at the top of each curve). For TEdd > 0.639 
the curve tends towards infinity when ui approaches w max 
(see Table 1). a is set to 0.43 (see text). 



M(fi) 
M(Q = 0) 



(1 



Edd J 



1 



iui 2 V 



Edd 



(16) 



To obtain Eq. (fTB")) . we used eq ((5]), and replaced the mean 
density p m by V/M, with V the volume enclosed by the 
stellar surface. 

For stars having an Eddington factor r E dd larger than 
0.639, we can use the expression of u cr it,2 in Eq. ^ to 
rewrite the last relation: 



M(ft) 
M(Q = 0) 



(1 



Edd J 



(i-r Edd ) i-- 



v 



(17) 



Thus, we see that the mass loss enhancement is governed 
by the ratio w/w C rit,2, while the deformation of the star is 
governed by the ratio w/w C rit,i- In the expression of u C rit,2 
intervenes the Eddington factor T^dd, thus the global en- 
hancement factor of the mass loss will depend on two pa- 
rameters UJ (or V C q/Vcrit,2) and I^dd- 

Fig. [5] shows the variation of the ratio M(ui)/M(uj = 0) 
as a function of these two parameters. For this plot, we 
took a value for a = 0.43, which is adapted for effective 
temperatures 4.05 < log(T c ff) < 4.3 (Lamcrs 2004, private 
communication) . The volume V enclosed by the stellar sur- 
face is numerically computed using eq. ((5J for the shape of 
the surface 

For Eddington factors greater than 0.639, w cr it,2 < ^crit,i 
(jMaeder fc Mevnetll2000l ). thus the second limit is the one 
to be considered. At this limit, the mass loss rate becomes 
very high. The precise value of the enhancement cannot 



Table 1. Maximum rotation parameter w max , and maxi- 
mum increase of the mass loss rate as a function of the 
Eddington factor I^dd (see text). 



TEdd 


^raax 


M(u)/M(0) 


TEdd 


^max 


M{uj)/M{Q) 


0.0 


1.0 


1.810 


0.6 


1.0 


21.696 


0.1 


1.0 


1.972 


0.639 


1.0 


oo 


0.2 


1.0 


2.214 


0.7 


0.968 


oo 


0.3 


1.0 


2.612 


0.8 


0.861 


oo 


0.4 


1.0 


3.383 


0.9 


0.659 


oo 


0.5 


1.0 


5.444 


0.95 


0.484 


oo 



be given since, at that limit and beyond, the hypothesis 
made in deriving eq. (fTB")) no long holds. For instance, near 
the Eddington limit, continuous radiation field contributes 
in pushing out the outer layers, while for obtaining the 
above expression, we made the hypothesis that the wind 
is triggered by radiation pressure on lines. 

In Table [JJ we indicate the maximum values of the rota- 
tion parameter w max above which the formulae given before 
for the enhancement of the mass loss rate due to rotation 
do no long hold. For values of T E dd inferior to 0.639, the 
maximum value is equal to one. In that range, stars with 
to reaching 1 will begin to lose mass through a mechanical 
mass loss in the equatorial regions. For T E dd superior to 
0.639, the maximum value is inferior to one. In that domain, 
the continuous emission will participate in pushing out the 
matter and very impo rtant mass loss rates are expected 
(jvan Marie et al.l 12008). Depending on the value of TEdd, 
the winds can be more or less anisotropic: for T^dd values 
just above 0.639, u cr it,2 is near f cr it,i and strong anisotropics 
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Fig. 9. HR diagram for the isotropic model (black solid 
line) and anisotropic model (red dashed line) with realis- 
tic mass loss rate. The rotation rate is indicated for some 
points along the tracks. 
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Fig. 10. Top panel: Rotation rate as a function of the cen- 
tral hydrogen mass fraction for the isotropic (black solid 
line) and anisotropic model (red dashed line) with realistic 
mass loss rate. The evolution proceeds from left to right. 
Bottom panel: Total angular momentum contained in the 
star in unit of 10 53 gcm 2 s _1 . The models are represented 
as in top panel. 



are expected; when I^dd is near 1, w cr it.2 is much lower than 
u c rit,i and the winds are expected to be isotropic. 

4.3. Models with realistic mass loss rate 

Here we discuss 9 models using a realistic mass loss 
rate (|Vink et al.ll200l[ ). and accounting for the increase of 
the mass loss rate induced by rotation (see above). One 
model takes the effect of anisotropic winds into account, 
and one has an isotropic mass loss. We also included in 
the anisotropic model the variation of the force multiplier 
parameters k and a due to the variation of the local effective 
temperature as a function of the co-latitude. Models were 
followed until they reached the critical velocity. 

The variations of the force multiplier parameters arc 
accounted for in the following way. In the expression of the 
increase of the global mass loss rate due to rotation ([TBI) , 
we take a and k given by the mean effective temperature 
of the star. For the computation of the anisotropic effects 
(JlJ, we use at each co- latitude a and k corresponding to 
the local effective temperature, allowing variations over the 
stellar surface. 

Fig. IH1 shows the HRD for this set of models. The black 
curve represents the model with an isotropic mass loss over 
the surface, and the red curve the anisotropic one. As in the 
previous case, the tracks in this diagram are very similar 
for both models, even if more complete physics is included. 
This confirms that even with full treatment of the wind 
anisotropy, including a realistic mass loss rate and the vari- 
ation of the force multiplier parameters, the effect of the 
anisotropic mass loss for very fast rotators remains very 
small. 



On the top panel of Fig. [TU1 we see the evolution of the 
ratio w as a function of the central hydrogen content for 
the isotropic model (black solid line) and the anisotropic 
one (red dashed line). The ZAMS is on the left, and the 
point where u> = 1 is on the right. The anisotropic model 
rotates ever faster than the isotropic one, but the difference 
is limited. The increase of rotation of the isotropic model 
near X c = 0.52 causes the crossing of the tracks in the 
HRD (see Fig. The inflation of the surface induced by 
the higher rotation rate decreases the mean effective tem- 
perature, pushing the track on the right. 

The bottom panel of Fig. [10] shows an interesting fea- 
ture. From the ZAMS until X c ~ 0.4, the behaviour of 
the total angular momentum contained in the whole star 
is similar to the models with constant mass loss rate: the 
anisotropic model keeps more angular momentum due to 
the wind being polar. From that point on, however, the star 
reaches a high enough rotation rate to strongly decrease 
the equatorial effective temperature. The force multiplier 
parameters are different in this area, and generate a strong 
equatorial mass loss. This effect produces the change in the 
angular momentum loss rate: the anisotropic model loses 
more angular momentum than the isotropic model. When 
the star reaches the critical limit, the isotropic model fi- 
nally has a higher angular momentum content than the 
anisotropic model! 

The final angular momentum of the star strongly de- 
pends on the angular momentum removed by the mechan- 
ical mass loss that the star undergoes during the criti- 
cally rotating phase. It is difficult to estimate which of the 
isotropic or anisotropic model will have the higher content 
at the end of the stellar evolution, and to quantify this 
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difference without a model accounting for the mechanical 
mass loss. This question, and first estimates of the mass lost 
in the equatorial disk, will be addressed in a forthcoming 
paper. 

5. Conclusion 

The main result of this paper is that radiative wind 
anisotropics do not strongly affect the angular momentum 
content of stars, in contrast with previous estimates. The 
different conclusion obtained here comes mainly from a pre- 
cise account of the effect of the surface deformation in ad- 
dition to the effects induced by the variation of the mass 
flux with the co-latitude. Interestingly, taking into account 
the variation of the force multiplier parameters over the 
surface when the star is near the critical limit can favour 
an equatorial-enhanced mass loss rather than a polar mass 
loss. In that case, the angular momentum loss when the 
effects of wind anisotropics are accounted for can be higher 
than when they are neglected! 

Since the anisotropic winds do have a strong influence 
on the evolution of the star, the strong enhancement of the 
polar mass flux has a b ig effect on the evol ution of the cir- 
cumstellar medium fsee lGeorgy et al1l2009l ). The formation 
of an asymmetric nebula around fast rotating stars is likely. 

However, another point which appeared in that work is 
the importance of the equatorial mass loss triggered by the 
reaching of the first critical limit. Two effects can act to 
keep the star at the critical limit: the first is the variation 
of the force multiplier parameters a and k in the equatorial 
regions when the local effective temperature becomes low 
enough due to the effect of rotation. This triggers strong 
equatorial radiative winds, and is already accounted for in 
this study. The other is the mechanical mass loss in the 
equatorial plane, when the equatorial effective gravity van- 
ishes. In that case, the mechanical mass loss through an 
equatorial disc can remove angular momentum. In a forth- 
coming paper we shall study in a quantitative way the im- 
pact of such a disc mass loss on the evolution of fast rotating 
stars. 
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